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Motivated by the recent experiments on Bose-Einstein mixtures with tunable interactions we study 
repulsive weakly interacting Bose mixtures at finite temperature. We obtain phase diagrams using Hartree- 
Fock theory which are directly applicable to experimentally trapped systems. Almost all features of the 
diagrams can be characterized using simple physical insights. Our work reveals two surprising effects 
which are dissimilar to a system at zero temperature. First of all, no pure phases exist, that is, at each point 
in the trap, particles of both species are always present. Second, even for very weak interspecies repulsion 
when full mixing is expected, condensate particles of both species may be present in a trap without them 
being mixed. 
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Introduction - Since the early realization of mixtures 
of Bose-Einstein condensates (BECs) ten years ago 
a host of experiments revealed the rich physics of two- 
component Bose 121, |3|] and Fermi systems. Whereas 
most experiments use mixtures of hyperfine states of hke 
atomic species, it is also possible for unlike species to be 
trapped simultaneously. Moreover both systems allow the 
use of Feshbach resonances by means of which the inter- 
atomic interactions can be tuned to arbitrary values and 
which revolutionized the physics of ultracold gases. Such 
resonances were observed for binary Bose gases consist- 
ing of 87Rb-4iK H, 87Rb-i33cs ii and ^^Rb-^^Rb fli 
mixtures; clear phase segregation was achieved for the last 
mixture by changing the intraspecies interactions of the 
s^Rb particles H. 

To date theoretical works on BEC mixtures predomi- 
nantly focussed on static and dynamic properties at zero 
temperature Interesting questions, however, can be 
posed concerning the finite temperature regime: for in- 
stance, will, upon increasing the interspecies repulsion, 
phase segregation for thermally depleted particles set in at 
the same time as for condensate particles? Does increasing 
the temperature induce an increased tendency to mix the 
species? Can pure phases still exist? The issue of phase 
segregation at nonzero temperature was already addressed 
in Ref. 10 where the authors calculated the conditions for 
an instability to occur in a volume with fixed particle num- 
bers. Also the case of a Bose mixtures in a trap at nonzero 
temperature was briefly discussed. 

In this Letter, we aim at clarifying the possible phase 
structures in traps of Bose mixtures at finite temperature. 
The presented phase diagrams are directly applicable to 



trapped situations, and yet, independent of the specific 
trapping parameters. We realize this by drawing the dia- 
grams as a function of the chemical potentials of the two 
species while taking the temperature to be fixed. These di- 
agrams are apt for gases confined by a smoothly varying 
potential for then the use of a local density approximation 
is appropriate and an effective chemical potential can be 
taken as locally constant at each position in the trap. Rather 
than giving exact diagrams using specific parameters [17], 
we focus here on their basic understanding by highlighting 
the generic features; more elaborate work will be presented 
in a forthcoming paper. We argue that our diagrams can 
be almost fully characterized, based on the knowledge of 
single-species Bose systems at finite temperature and the 
BEC mixture at zero temperature. 

Equation of state — Consider two Bose species, la- 
belled 1 and 2, at temperature T and at chemical poten- 
tials fi^ and The particles of species i and j inter- 
act weakly via s-wave scattering, quantified by a posi- 
tive scattering length a., and a coupling constant G.. = 
2TTfi^a..{m1[^ + mj^) (henceforth i,j = 1,2). We use 
the self-consistent Hartree-Fock (HF ) equations of state as 
first introduced by Huang et al. lllin . This very model is 
known to give good agreement with experiments and 
even with quantum Monte Carlo simulations [13J. The HF 
model treats noncondensed particles as free particles with a 
mean field chemical potential shift, and takes only thermal 
and no quantum fluctuations into account. The HF grand 
potential per unit volume Q is expressed in terms of the 
condensate densities n^. and thermally depleted densities 



n 



E 

i=l, 2 



^5/2 



A3 



n 



2 



where A, = j (m.kgT), P = l/k^T and g,{x) = 

s^'x'. The first two terms are due to the entropy and 
the kinetic energy of the depleted particles. The very last 
term couples the two species by repulsive interactions. At 
fixed chemical potential, minimization of 17 with respect to 
the densities n^. and n^. yields the self-consistent Hartree- 
Fock equations of state: 
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and, 11° = —G..n^. whenever n^^ ^ 0. The depleted 
densities on the other hand are governed by minimization 
with respect to IJ,° : 
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For most systems of ultracold gases only few length- 
scales are relevant. These include the particle de Broglie 
wavelength A., the interparticle distance n"^/'^ and lastly 
the scattering lengths a... The HF theory as introduced 
above is valid in case of weak interactions i.e. for small 
values of a.. / A . and when ^/ng/e -particle excitations pre- 
vail i.e. when n.a..X^ <C 1 [2]. Also, the theory fails 
very close to the transition point where fluctuations become 
dominant. In nowadays experiments, the scattering length 
takes values of order of ten nanometers and below, whereas 
the interparticle distances and the Broglie wavelengths are 
typically of the order of microns and higher. It follows that 
the necessary conditions for the HF theory to be valid are 
fulfilled m. 

We proceed by first discussing BEC mixtures at zero 
temperature and a single species Bose gas at finite temper- 
ature, since these systems provide all ingredients necessary 
for the understanding of the results obtained further. 

BEC mixture at T = — At zero temperature, the HF 
theory reduces to Gross-Pitaevskii (GP) theory for a spa- 
tially homogeneous system. According to GP theory, two 
regimes can be distinguished, depending on the parameter 
A defined as H: 
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When A < 0, the inter-species repulsion is sufficiently 
strong to induces phase segregation and solely BEC 1 or 
BEC 2 may exist as ground states, that is, species 1 and 2 
are immiscible. One readily finds the following conditions: 



BEC 1 when > /i, ^fCJCT^ 



BEC 2 when /^^ < /i, ^JGJG^ 
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On the other hand, when A > 0, the interspecies repulsion 
is weak and a mixed phase of nonzero density of BEC 1 
and BEC 2, here denoted as BEC MIX, may appear as the 
ground state: 

BEC MIX when G,, /G,, < /x^ < G^JG^^ . (5) 




FIG. 1: (Color online) Generic phase diagram for two Bose 
species at fixed temperature in case of strong interspecies repul- 
sion (A < 0) as a function of the chemical potentials of both 
species. The BEC 1 (green) phase consists of condensate parti- 
cles of species 1 and thermally depleted particles of species 1 and 
2 and vice versa for the BEC 2 (blue) phase. The (grey) phase NO 
BEC consists of depleted particles of both species. The chemical 
potential /if is given by Eq. Q and the locus of the triple point 
TP is given in Eq. (|7]i while the first-order phase boundaries are 
described by Eqs. ^ and dlOl l. 



Single-species at finite T — Consider now a gas of 
species i at fixed temperature. One can think of the chem- 
ical potential /i. as a measure for the particle density: the 
higher /i., the higher the density of species i. For large 
and negative values of /i . , the system is dilute and the in- 
teraction energy small. On the other hand, above a cer- 
tain "critical" positive chemical potential a transition to a 
phase with nonzero condensate density sets in; by the HF 
approach, this critical value is attained when = such 
that H CI: 

= 2G,, C(3/2)/A3, (6) 
with C the Riemann zeta function. Note that nf de- 
pends on temperature via A. and from the scaling /if ~ 
kgT{a.J\.) we conclude that /x^ <^ k^T. 

The introduction of a second Bose gas will shift the criti- 
cal chemical potential to a higher value than /i^ ; this can be 
ascribed to the interspecies repulsion which tends to rarify 
the gas of species i. 

Phase Diagram for A < and finite T — Upon work- 
ing at fixed particle numbers as in Ref. 10, the onset of 
phase segregation can be naturally probed by the presence 
of a density instability. At fixed chemical potentials on the 
other hand, this method is inadequate as an instability in- 
dicates a spinodal (line) rather than the transition itself in 
case it is of first order. Instead, we have compared numer- 
ically the grand potentials for the eight possible thermody- 
namical phases and their associated multiple solutions. We 
arrive at the generic phase diagrams for Bose mixtures at 
finite temperature in Figs. [T] and |2] for the cases A < and 
A > respectively. 
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Identical to the behavior at T = 0, no mixing of unlike 
condensate particles is possible at finite T when A < 0. 
Only three phases appear in Fig.[T] none of which are pure; 
that is, even for extremely strong interspecies repulsion, 
depleted particles of both species are present in all phases. 
Accordingly, BEC 1 denotes the phase composed of con- 
densate particles of species 1 and depleted particles of both 
species, and vice versa for BEC 2. The NO BEC phase de- 
notes the phase with thermal particles of species 1 and 2. 

A triple point TP is present in the phase diagram of 
Fig. [TJ As checked numerically, its locus (/i^^ , /i^^ ) can 
be found by setting fj,'^ = ^u" = in Eq. ([T]l: 



llY = ^l- + ^x<;G,J2G,,. 
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Note that in units of /i^ and //J, the position of point TP is 
temperature independent. 

Consider now the phase boundary between BEC 1 and 
2; for large and positive values of both chemical potentials, 
the interaction energy per particle will eventually become 
much larger than the thermal energy and so asymptotically 
the T = behavior is expected. Indeed, we find that in 
accord with Eq. (|4]l the BEC 1-2 phase boundary satisfies: 



^1 //^2 = yGJG^^ when both [i^ , 
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Moving along the NO BEC-BEC 1 boundary away from 
TP, the density of species 2 decreases; in the limit of 
fi^ — > — oo, only particles of species 1 remain such that the 
single-component value fi^ fi'^ is attained. The associ- 
ated asymptotic decay is exponential as the NO BEC-BEC 
1 phase boundary is described by: 
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Because ^u^ <C T, it follows that the asymptotic conver- 
gence is slow if ^2 is expressed in units of /i^ . The same 
line of argument with interchanged indices 1 and 2 applies 
to the demarcation line between BEC 2 and NO BEC: 



/uf + G,,e^i/'^B^/A^ when fi^ -oo. (10) 



In outline, we have shown that one can qualitatively es- 
tablish the phase diagram for the case A < without ex- 
plicitly performing any numerics. Indeed, we found that 
the expressions for the triple point (Eq. (|7]l), and the phase 
boundaries (Eqs. (l8]|-([Toll) agree well with the numerically 
exact results. 

Phase Diagram for A > at finite T — At zero tem- 
perature mixing of unlike condensate particles is possible 
when A is positive. We find the same behavior at finite 
T. This is shown in Fig. |2] where BEC MIX consists of 
condensed and depleted particles of both species. The bi- 
furcation of the BEC 1-2 line brings along the appearance 
of a new triple point TP2 (TPl is identical to TP of Fig.[T]l. 
For the aforementioned reasons, the T = physics is re- 
covered for large and positive chemical potentials. Hence 
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FIG. 2: (Color online) The same applies as in Fig. [T] except that 
weak interspecies interactions are considered (A > 0). The (dark 
grey) BEC MIX phase consists of condensed and depleted par- 
ticles of both species. The red arrowed line indicates a possible 
tracing of the diagram as a function of the (radial) position in the 
trap. The point B indicates the chemical potentials at the center 
of the trap. 



Eq. (|5]l implies that the BEC MIX phase is delimited by the 
phase boundaries 



^l//^2 = ^^12/^*22 



when both , ^ 00. (11) 



The location of triple point TP2 alters upon varying tem- 
perature and A. For small and positive values of A, TP2 
enters at large and positive values of the chemical po- 
tentials while TP2 coalesces with TPl only in the limit 
A 00. This leads us to a remarkable conclusion: even 
when the zero temperature condition for species mixing is 
satisfied, at finite temperature, two-phase equilibrium be- 
tween BEC 1 and 2 is still allowed. Changing the temper- 
ature also induces a displacement of TP2 with respect to 
TPl. Intuitively one expects that increasing the tempera- 
ture favors mixing. Yet, condensate particles carry no en- 
tropy and we find rather the opposite: increasing the tem- 
perature suppresses mixing. In particular, the distance be- 
tween TPl and TP2 (in units of /i^) increases upon increas- 
ing the temperature. 

In the ideal case of fully suppressed interspecies inter- 
actions (A — > 00) the triple points TPl and TP2 coincide 
and, as expected, the phase boundaries are the horizontal 
and vertical lines /i^ = //^ and =1^1- 

To summarize, the A > phase diagram is similar to 
the A < diagrams but is marked by the appearance of an 
additional phase consisting of mixed BEC species, and an 
additional triple point TP2. Despite the known asymptotic 
behavior of the phase boundaries around the BEC MIX 
phase (see Eq. (fTol)). the locus of TP2 depends in a non- 
trivial way on temperature and A and must therefore be 
determined numerically. 
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Discussion — In order to extract from Figs. [T] and |2] 
the possible phase structures appearing in an experimental 
trap, the trapping potentials U. are required. In most ex- 
periments a local density approximation or Thomas-Fermi 
approximation is justified such that at each position r an 
effective chemical potential /i. (r) = fi. — U. (r) may be 
assumed constant; /i . (r) is maximal at the center and min- 
imal at the edge of a trap. In case both species are har- 
monically confined or U.{r) = m.uP'.v^ l2, the chemical 
potentials probe our phase diagrams following a straight 
line: 

^^^lii^A"^) - fJ-i) = m,u!^^{ii,ir) - /ij. (12) 

As an example, a trap configuration which is described by 
the red path in Fig. |2] contains a BEC 2 core surrounded 
by shells of BEC 1 and NO BEC. Note that our phase dia- 
grams do not explicitly depend on trapping parameters nor 
gravity. 

The most striking conclusion of this work is that no 
pure phases exist at finite temperature; even for large inter- 
species interactions, the depleted particles are not entirely 
expelled |20]. However, the fact that upon increasing the 
interspecies repulsion G^^ in the BEC MIX phase, con- 
densate particles will be expelled first, may be understood 
from the following simplistic argument: in the expression 
for il, it is seen that due to different exchange terms the in- 
traspecies interactions for depleted particles are a factor 2 
higher than those of condensate particles. A naive applica- 
tion of the zero-temperature criterion for phase segregation 
> G^^G^^) to depleted particles, yields segrega- 
tion for values of G^^ higher than ^JIG^^G^^. 

As a direct extension of this work, the interface physics 
and its impact on the phase structures can be explored in 
case of two-phase equilibrium. A generalization of our HE 
theory to spatially inhomogeneous systems is thereby re- 
quired; within a first approximation, this can be effectu- 
ated by introducing the terms V^y^rT" into the grand po- 
tential [3]. Such theory would then allow the calculation of 
surface excess quantities which are useful for tightly con- 
fined gases where their effect on the ground state configu- 
ration turns out to be substantial ll8i[l4l[l5ll. Further chal- 
lenges include a generalization of the theory of anomalous 
wetting phase transitions to finite temperatures [16]. 

Conclusion — We have elaborated generic phase dia- 
grams for Bose mixtures at finite temperature based on the 
Hartree-Fock model. As corroborated by numerical analy- 
sis, almost all properties of our phase diagrams can be ex- 
pressed in terms of known results for single-species gases 
at finite temperature and for binary mixtures at zero tem- 
perature. This means that it is possible to establish the 
phase diagrams (aside from the triple point TP2 in case 
A > 0) without performing any numerics. The phase dia- 
grams are drawn as a function of the chemical potentials of 



the species such that, spatially tracing the trap is equivalent 
to exploring the phase diagram. We highlight the impor- 
tance of the triple points in our diagrams and find that, sur- 
prisingly, no pure phases exist, and, increasing the temper- 
ature tends to suppress the mixing of condensate particles. 
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